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Motivated by the recent discoveries of spin-1 and spin-1/2 
Bose gas, we have studied the general structure of the Bose 
gases with arbitrary spin. A general method is developed to 
uncover the elementary building blocks of the angular mo- 
mentum eigenstates, as well as the relations (or interactions) 
between them. Applications of this method to Bose gas with 
integer spins (/ — 1, 2, 3) and half integer spins (/ = 1/2, 3/2) 
reveal many surprising structures. 

Recent experiments on dilute quantum gases of alkali 
atoms have produced a spin-1 and (pseudo-) spin-1/2 
Bose gas respectively. The former is produced in opti- 
cally trapped 23 Na jjj, the latter in magnetically trapped 
87 Rb by rotating the hyperfine states \ f = 2, m = 1) and 
|/ = l,rn = —1) into each other through a slightly de- 
tuned rf-field Q. Many macroscopic quantum phenom- 
ena have been observed in these systems. At present, 
these phenomena can be explained in terms of the sin- 
gle condensate picture. However, in the case of spin-1 
Bose gas with antiferromagnetic interaction like 23 Na, 
it is been pointed out very recently that as the magnetic 
field gradient is reduced, the single condensate will evolve 
toward an angular momentum eigenstate, which will be- 
come a spin singlet as the magnetic field is reduced to 
zero H M . The singlet state is a "fragmented" structure 
which bears no resemblance to the single condensate state 
j| . That the ground state of a Bose system can be very 
different from a conventional single condensate when it 
acquires internal degrees of freedom is a surprise. 

Motivated by the fragmented structure of the spin-1 
Bose gas, we consider Bose gases with higher spins. Al- 
though Bose gases with spin / > 1 have not yet been 
produced, it is conceivable that they can be realized in 
the future. After all, both spin-1 and spin-1/2 Bose gases 
have only come into existence within the last one and 
half year. Currently, there are also efforts to condense 
85 Rb, which will be a spin-2 Bose gas when loaded into 
an optical trap. The main reasons for our investigation, 
however, remains theoretical and conceptual. The nature 
of the ground states of Bose gases with internal degrees 
of freedom is of fundamental importance. It has a place 
in the lore of superfluid physics and significance that goes 
beyond to the study of Bose-Einstein condensation. Our 
goal is to present a general method to construct the (to- 
tal) angular momentum eigenstates \F,F Z ) for Bose gases 
with arbitrary spin /. The construction of these eigen- 
states is a crucial step in diagonalizing the Hamiltonian of 
the system. Our method reveals many surprising struc- 
tures. Generally, the spin state \F, F z = F) is made up 
of singlet and magnetic building units. A schematic rep- 
resentation of the structure of the spin state \F, F) for 
Bosons with spins / = 1,2,3 and "pseudo-spin" 3/2 are 



shown in figure a to d. They illustrate the intricate struc- 
ture of these eigenstates and their increasing complexity 
with increasing /. 

The essence of the problem can be illustrated by con- 
sidering a homogeneous Bose gas with spin-/. (Its rela- 
tion to a trapped gas can be understood cither in terms 
of local density approximation and in the procedure out- 
lined in ref. Q.) For a homogeneous dilute Bose gas, 
we first consider the condensate in the zero momentum 
mode (i.e. k = 0), denoted by the annihilation operator 

= a M (k = 0), where \x labels the 2/ + 1 spin compo- 
nents. The angular momentum operator then becomes 
F = ajjf^et^, where f Mt/ is the spin matrix for a spin-/ 
Boson. The effect of the k ^ modes is to deplete the 
condensate. However, they can be ignored in the zeroth 
order approximation as they only contribute a small cor- 
rection to the energy. (For trapped gases, the k = 
mode will be replaced by the lowest self-consistent mode 
that the system condenses into || Q). To construct the 
angular momentum eigenstates, it is sufficient to focus on 
the states \F, F z = F) with maximum spin projections, 
since other states with F z < F can be obtained by apply- 
ing to \F, F) the spin lowering operator F_ = F x — iF y . 
In the following, we shall first derive our method, and 
then illustrate its application for the integer cases / = 1 
to 3 and half integer case / = 1/2 and 3/2. The case 
/ = 3 is particularly subtle and will be considered last. 

(LI) Outline of the Generating Function Method: To 
give an orientation of our discussions, we first outline 
the logic of our method before presenting the detailed 
derivations. We begin by considering the total number of 
maximum spin states \F, F) for a systems of N particles, 
which we denote as Mn(F). To generate this number for 
all N and S simultaneously, we consider the generating 
function 



G(x,y)=J2T, M * 



{F)x N y F . 



(1) 



N>0 F>0 



where x and y are complex numbers within the unit circle 
(|x|, \y\ < 1) to ensure convergence. Once this function 
is constructed, we shall see that Mn(F) is given by the 
number of solutions of a set of equations obeyed by two 
sets of non-negative integers {si > 0} and {rrij > 0}. The 
integer Si is the number of singlet building unit O^ which 

(s) 

is made up of n\ Bosons and carries no spin, while mj is 
the number of magnetic building unit Tj which is made 

up of Bosons and carries spin £j. The integers {s; > 
0} and {mj > 0} satisfy number and spin constraints 



E 



(m) 

Uj rrij 



= F. 



(2) 



1 



as well as a set of conditions C a that further limit the 
range of the integers {s^ and the {rrij}. These condi- 
tions C a reflect the inter-dependence (or "interactions" ) 
among the building units. The conditions C a are very 
simple for spin / < 3 but become quite complicated as 
/ > 3, illustrating the rapidly increasing complexity of 
the system as / increases. The typical form of these 
conditions will become clear when we come to our exam- 
ples. The general structure of the maximum spin state is 
therefore \F,F) = £ A({ Si }, {rrij}) U itj eJ*Tj m ' |vac>, 
where the A's are coefficients and the sum is over all 
non-negative integers {s,; > 0} and {rrij > 0} satisfying 
the constraints C a . 

(1.2) Derivation of the Generating Function Method : 
We begin with the observation that the integer Mn(F) 
can be expressed as 



M N (F)=I N (F)-I N (F + 1), 



(3) 



where In(F) is the total number of states with F z = F, 
independent of the value of total spin F. Eq.(|^) follows 
from the fact that all spin multiplets with total spin F' > 
F will contain a state \F',F Z = F), which contribute 1 
to both I N (F) and I N (F + 1), and hence to M N (F). 
Only those spin states with total spin F total = F z = F 
will be included in the integer In(F) and not In(F + 
1). That eq.(|) is useful is because it is much easier 
to construct a generating function for In(F) due to the 
removal of the spin constraint. Before proceeding, we 
note that while M N (F) is defined only for F > 0, In(F) 
is defined for both positive and negative F such that 
I N (F) = I N (-F). 

To find In(F), we note that a many-body state with 
total spin projection F z = F is of the form 



\F,F, 



F) 



with y% = . 



N 



m {,,,})[ J] ap\ |vac) (4) 
F, where {nf\ is 



E 

{nj>0} \j=-f 



-f 



-. <=-f n 3 = iv and Hi - ' vv^c^ v'jj 

a set of 2f + 1 non-negative integers, a] creates a Boson 
in spin state j, and -B's are coefficients. The number of 
states with F z = F is 



IN 



(F)= E 



{nj>0} 




F 



(5) 



where A(x) is a delta- function ensuring the vanishing of 
x. The generating function eq.(Q) can in principle be ob- 
tained by substituting eqs.(||) and ^ into eq.(|l]). How- 
ever, the constraint F > in eq.(|l]) prevents an efficient 
summation. We therefore consider the function 



w(x, y) = EE - In ( f + !)) xN y F > 



(6) 



A'>0 F 



where the sum F ranges over all integers. Clearly, G(x, y) 
is W{x, y) with all negative powers of y eliminated. This 
elimination can be achieved by the following integration 



r27T in 00 

G (^y)= / ^E[(^~ l ) v (^)] 

Jo 1=0 



(7) 



Performing the sum in eq. (^) , G becomes a contour inte- 
gral around the unit circle C, z = e lB , 



G{x,y) 



dz W(x, z) 
c 2iri z — y 



(8) 



The expression W can be obtained easily since F now 
runs through all integers. Substituting eq.(||) into eq.(||), 
and first sum over F and N, the functions W becomes 



i=-f 




(9) 



(10) 



We then arrive at the key expression for the generating 
function, 



G(x,y) 



r dz 1 - z- 1 -pj f 
J c 2iri z-y .^1 1 - x& 



(11) 



To illustrate how eq. (jlfj) can be used to obtain the struc- 
ture of the maximum spin state \F,F), we consider the 
following examples. 

Spin-1 Bosons : For / = 1, eq.(O) gives 



1 



{l-x 2 ){l-xy) 

= E E^+V 1 . 

n 2 >0^i>0 

Comparing with eq. (Q) , we have 

m n (s) - E E A ( 2nz + li - N ) A ^ 

n 2 >0^i>0 



(12) 
(13) 

F) ■ (14) 



Eq.(^4[) shows that Mn(S) is the number of the solutions 
of the equations 



2 • n 2 + 1 • h = N, 1 • 4 = F. 



(15) 



It is clear that eq. ([if]) has a unique solution no = (iV — 
F)/2,£ 1 = F. Hence M { J =1) (F) = 1, i.e. there is only 
one maximum spin state \F,F Z = F). Next, we recall 
that the exponent of x and y are associated with particle 
number and spin respectively. Eq.(|l5|) shows that the 
system consists of l\ magnetic structural units which are 
spin-1 Bosons (ai), and n 2 singlet pairs of Bosons 9 2 . A 



simple exercise shows that 8 2 = (2aia_ 



The (un- 



normalized) many-body state \F, F z = F) is then given 
by \F,F) = a\ F Q^ N ~ F ^ 2 \\ac) , which is the results given 
in ref.( j§). (See also fig.a). 
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Spin-2 Bosons : For / = 2, eq.(ll) gives 



G^= 2 \x,y) 



3, ,3 



1 + x A y 



(1 - x 2 ){\ - x 3 )(l - xy 2 ){l - x 2 y 2 ) ' 



(16) 



'y ' ^ ' ^2s2+3s3+mi+2m2+3m3y2mi+2m2+3m3 (17) 



m 3 =0,l 



where the first sum is over the non-negative integers set 
{sg, S3, mi, ma}. We then see that Mn(F) is given by 
the number of solutions to the equations 



2s 2 + 3s 3 + mi + 2m 2 + 3m 3 = N, 



2mi + 2m2 + 3rn3 = F. 



(18) 



(19) 



A solution of eqs.(|l8|) and ( |l9| ) describes a state consist- 
ing s 2 two-particle singlets 82, S3 three-particle singlets 
03, mi spin-2 Bosons (a 2 ), and m2 two-particle spin-2 
state 1 2, 2) (denoted as T 2 ). Since m 3 = and 1, the sys- 
tem may or may not contain a three-particle spin-3 state 
1 3, 3) (dented as T3) depending on whether F is odd or 
even. It is straightforward to work out the expressions 
of these states, which are O2 = &2Q-2 — <Xi<X-i 

6 3 = a (2a 2 a_ 2 + a x a_ x - |a§) - 



ana, 



2"0 



and r.q 



2a2a_ x 



In 2 

2 a 0> 

2 (afa„ 2 + a 2 a 2 _ 1 ), 
\f&a 2 a l a a + af . 



The general structure of the state \F, F) is then \F,F) — 

EA({si},{m}) (a^r^r^) e 2 S2 e^|vac). 

The condition m3 = 0, 1 is the additional constraint 
C a mentioned in section (1.1). If r 3 was a "free" unit 
that could appear as many as times as possible, the nu- 
merator of eq.(^6|) would be (instead of 1 + x 3 y 3 ) an infi- 
nite series Y^m 3 =o( x3 y 3 ) m3 > which will turn into a factor 
(1 — a; 3 ?/ 3 ) -1 like other "free" building units (62,83,02, 
and F 2 ) in the denominators in eq.(jlq). The fact that 
the series of x 3 y 3 terminates at the first order means 
that a pair of three-particle singlets can be expressed in 
terms of all other "free" excitations (0 2 , 3 , a 2 , and r 2 ) 
and therefore has already been accounted for in the gen- 
erating function. Indeed, when examining T 3 (because 
of the prediction of the generating function), one finds 

r 2 = _ 16^6 p3 + 8^2 a 2p q 4 /| a 3 03 |||_ Not6i how _ 



3 — ~ x 2 "T" — ~ "2 X 2^2 — 3 «2 W 3 

ever, that I^ appears at most once. It therefore has no 
thermodynamic significance. This means that one can 
obtain the relevant thermodynamic structure by taking 
any term in the numerator of eq.([l6]). (See also fig.b). 

(II). Bose gas with half integer spins : When / is a half 
integer, it is useful to consider the generating function 



G(x,y) 



EE* 

JV>0F>0 



(F)x N y* F . 



(20) 



Proceeding as the integer case, the function W{x, z) in 
eq.@ becomes W(x, z) = 0—z- 2 )Uj=-fO--xz 2j )- 1 . 



Since / is a half-integer, W consists of even and odd 
powers of z. Using the previous method to project out 
all the negative powers in y, we have 



G(x,y) 



dz 1 — z 



r 2iti z — y - 1 A 1 — xz 2 ^ 
° y j=-f 



f 1 

n.T^ 



(21) 



Spin 1/2 Bosons : For / = 1/2, we have 

G (/=1/2) (x,y) = - J— = 5>*V) n/2 (22) 
1 — xy z — ' 



;;>0 



Since the equation n = N, n/2 = F has only one solution 
and forces F — N/2. This means that the total spin of a 
systems of spin- 1/2 Bosons is fixed by the particle num- 



ber N to be F = N/2, and \F,F = N/2) = </ 2 |vac). 
The system can be referred to as a "statistical ferromag- 
net" since the ferromagnetism is forced by statistics. 
Spin 3/2 Bosons : For / = 3/2, we have 



G^ 3 ' 2 \x,y) = 



1 + x A y 



3, ,3 



(l-x 4 )(l-xy 3 )(l-x 2 y 2 ) 



E E 



4s4+mi+2m2+3m3 



(y 2 ) 



(23) 



. (24) 



m 3 =0,l 



where the first sum is over non-negative integers 
S4,mi,TO2. The number of spin state \F,F) is given by 
the number of the solution of 

3 3 
4s4 + mi + 2mi + 3m,3 = N, 2 TOl 7712 ~2 n3 = ^' 

(25) 

which describes a state consisting of S4 four-particle sin- 
glets 64, mi spin-3/2 Bosons (i.e. a 3 / 2 ), and m2 spin-1 
pairs |1, 1) made up of two spin-3/2 particles (denoted as 
r 2 ). Since 7713 = and 1, the system may also contain 
a spin-3/2 three-particle state ||,|), (denoted as r 3 / 2 ) 
which appears at most once. Thus, we have \F, F) = 

£A({ Si }, W) (4?2 r 2 m2 r 3 m3 ) ei S4 |vac). (See also 

fig.c). 

(Ill) Spin-3 Bosons : The case of / = 3 begins to 
illustrate the full complexity of the Bosons with higher 
spin. It is sufficiently intricate so we discuss it last. When 
/ = 3, eq.(§) gives 

1 ' V> ~ (1 - x 2 )(l - x 4 )(l - a;6)(l - x™) [ ' 



D{x,y) 



(1 - xy 3 ){\ - x 2 y 2 ){\ - x 2 y*) 



(27) 



The term C{x, y) is a polynomial with about fifty terms 
of the form x a y b with (a, b > 0). Since b > 0, these terms 
represent magnetic structures. From eq. ([!]), we see that 



3 



the structure of the total singlet state |F = 0,F Z = 0) 
is given by G(x,y = 0). Extracting G(x,y = 0) from 
eq.(|l^) (i.e. setting C = and D = 1), we see that 
the singlet state is a linear combinations of singlets con- 
sisting of two, four, six, and ten particles, denoted as 
02, 64, 06, © 10 respectively. From our discussions for the 
spin-2 case, we see that all singlet except that made up of 
15 particles (©15) can be expressed as products and sums 
of the "free" singlet set {0 2 , ©4, ©g, ©io}- However, two 
15-particle singlets (i.e. @f 5 ) is reducible to free singlet 
units. (See also fig.d). 

As before, the elementary magnetic units {I\} are 
given by the denominator of D. They are single particle 
spin-3 Bosons (a 3 ), two-particle spin-2 pairs (|2, 2)), and 
two-particle spin-4 pairs (|4, 4)). The major difference 
between the / = 3 and previous examples, however, is the 
appearance of large number of terms in the numerator of 
the generating function (i.e. C), and the fact that about 
half of these terms have negative signs, which means dis- 
appearance rather than appearance of a configuration. 
The origin of the negative terms is due to the fact that a 
product of two or more different magnetic units T t and 
Tj can be expressed in terms of other magnetic and non- 
magnetic units. These are the "interaction" constraints 
C a we mentioned in section (LI), Note that in the case 
of / = 2, the interaction constraints comes from the re- 
ducibility of a single type of structure, i.e. T§ is reducible 
into other free units. As a result, all terms in the numera- 
tor of G( /=2 ) are positive because one simply enumerates 
the multiplicity of T 3 until it becomes reducible. If, how- 
ever, the interaction constraints involve the reducibility 
of the products of two or more different magnetic oper- 
ators, as well as "scattering" such as TiYj —t T^T/j+ctc, 
then the counting process can be not be simply a ter- 
mination of the multiplicity of a particular pattern. We 
shall not analyse the interaction constraint for the / = 3 
case here because it is very involved. Despite this com- 
plexity, it is clear from the generating function what the 
elementary magnetic building units are. 

In summary, we have illustrated the method to uncover 
the elementary building units of the angular momentum 
eigenstates of a spin-carrying Bose gas. The construc- 
tion of these units is crucial for energy studies JfJ. Yet 
even without such studies, the present method has illus- 
trated the complex structure of the ground state of these 
Bose gases as a function of magnetization. The fact that 
the number of independent singlet units proliferates as / 
increases also means that the system becomes more frag- 
mented, since spin fluctuations (which is already huge in 
the spin-1 case in the low field limit will increase as 
the number of different singlets increases. 
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Figure Captions : Figures a,b, and c are schematic rep- 
resentations of the basis of the angular momentum state 
\F, F) for spin / = 1,2, and 3/2 respectively. Enclosed 
units with and without arrows represent magnetic and 
singlet units respectively. The number of dots indicates 
the number of particles in the unit. For example, in the 
spin-2 case (fig.b), the state consists of 2-particle and 3- 
particle singlets (represented as arrow-free ellipses and 
triangles containing two and three dots resp.), and a 2- 
particle spin-2 pair (represented as an ellipse with two 
dots and an arrow). The 3-particle spin-3 unit is repre- 
sented as a triangle containing 3 dots and an arrow. The 
dashed circles in the interior are drawn to help to visu- 
alize the singlet and the magnetic units. They are not 
meant to imply the existence of a singlet core. Figure 
d is a schematic representation of the singlet structure 
of a spin-3 Bose gas, which consists of 2-, 4-, 6- , and 
10-particle singlets, and a "constraint" unit consisting of 
15 particles, which is reducible to other existing singlet 
units when appears more than once. 
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